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Abstract — The problem of open-loop dynamical control of 
generic open quantum systems is addressed. In particular, I 
focus on the task of effectively switching off environmental 
couplings responsible for unwanted decoherence and dissi- 
pation effects. After revisiting the standard framework for 
dynamical decoupling via deterministic controls, I describe 
a different approach whereby the controller intentionally 
acquires a random component. An explicit error bound on 
worst-case performance of stochastic decoupling is presented. 

I. Introduction 

The need for accurately controlling the dynamics of a 
quantum-mechanical system is central to a variety of tasks 
ranging across contemporary physics, engineering, and in- 
formation sciences [1], [2], [3]. In particular, motivated 
by both continuous experimental advances in nanoscale 
devices and the challenge to practically implement fault- 
tolerant quantum information processing, control strategies 
for open quantum systems undergoing realistic irreversible 
dynamics [4] play an increasingly prominent role. 

Dynamical decoupling techniques offer a versatile control 
toolbox for open quantum-system engineering [5], [6], [7]. 
In its essence, a decoupling protocol consists in a sequence 
of open-loop transformations on the target system (control 
pulses in the simplest setting), designed in such a way that 
the effect of unwanted dynamics is coherently averaged out 
in the resulting controlled evolution. Applied to the removal 
of unwanted couplings between the target system and its 
surrounding environment, this paves the way to a general 
strategy for decoherence control and error-suppressed quan- 
tum computation purely based on unitary control means. 

Both within formulations of the decoupling problem and 
more general coherent-control settings, the restriction to 
purely deterministic control fields has provided a most 
natural starting point. In a way, this finds ample justification 
in the fact that non-deterministic effects (such as stochastic 
noise and/or random control imperfections) typically dete- 
riorate system performance, motivating the effort for de- 
signing intrinsically robust decoupling schemes [8] and for 
assessing open-loop fault-tolerance thresholds [9]. Yet, no 
fundamental reasons exist for not lifting such a restriction, 
by purposefully allowing stochasticity in the underlying 
control design. Beside being conceptually intriguing on its 
own, it is worth recalling that notable examples may be 



found of situations where noise and randomness might have 
a beneficial rather than detrimental effect. Of special rele- 
vance are phenomena like the self-averaging of intermolecu- 
lar interactions in gases and liquids via random microscopic 
motions [10] and quantum stochastic resonance [11], or the 
idea of dissipation-assisted quantum computation [12]. 

A first step toward exploring randomized quantum control 
was recently taken by Viola and Knill [13], confirming in 
principle the possibility of enhanced system performance 
as compared to deterministic control in relevant scenarios. 
It is the purpose of this paper to further elucidate the 
random decoupling framework, by first presenting a general 
control-theoretic formulation and contrast it to the standard 
deterministic one (Section II), and then discuss in detail a 
quantitative error bound on stochastic control performance 
(Section III). Final remarks conclude in Section IV. 

II. Formulation of the control problem 

A. Quantum- control systems 

The standard open-loop control problem for an isolated, 
closed quantum system S defined on a state space Tis of 
dimension ^5 < oo is described (in units where ?i = 1) by 
a bilinear control system of the form [14] 

in 

H,(t) = Y,Hiu,(t). (1) 

1=1 

Here, U{t) is the evolution operator (or propagator) of the 
system, whereas _ffo = Hs, Hg represent the internal (or 
drift) Hamiltonian, and the applied control Hamiltonians, 
respectively. Both Hq and the Hi are Hermitian operators 
on Jis which, without loss of generality, may be assumed 
to be traceless. The time dependence of the overall control 
Hamiltonian Hc{t) is modeled through the real functions 
ui{t), which typically represent electromagnetic fields and 
are the control inputs of the problem. A broad separation 
between deterministic and stochastic control systems may 
be drawn depending on whether each control input is 
a deterministic function of time or some randomness is 
allowed for at least one input. The state of S is described 
in general by a Hermitian, positive operator ps on Tis, 
normalized with respect to the trace norm in such a way 



that trs{ps) = 1- In what follows, I will assume that S 
is initially in a pure state, described by a one-dimensional 
projector tts of the form tts — \il>){'ip\, with £ Hs- 

It is convenient to focus directly on the control propaga- 
tor Uc{t) as the basic object for control design, 



Uc{t) = r exp 



duHc{u] 



(2) 



where the symbol T denotes as usual time ordering. By 
effecting a canonical transformation to a time-dependent 
frame that continuously follows the applied control. 



psit) = ui{t)psms) , 



(3) 



the explicit action of the control field is removed from the 
dynamics. The control problem of Eq. Q takes the form 



dU{t) 
dt 

Hit) 



-iH{t)U{t) , 
Ul{t)HoU,{t) 



(4) 



in terms of the propagator lJ{t) for the transformed state, 

Ps[t)^U{t)ps{mHt) . U{t) = Ul{t)U{t) . (5) 

I will refer to the formulations of Eqs. Q, (0} as physical 
and logical frame formulations, respectively. While from the 
mathematical point of view the logical frame description has 
the disadvantage of being highly non-linear in the control 
inputs, Eq. (|4} makes it very convenient to directly map 
properties of the desired effective evolution back into design 
constraints for Uc{t), and viceversa. If the control strategy 
is cyclic, that is Uc{t + Tc) = Uc{t) for Tc > 0, and Hq is 
time-independent as assumed so far, the periodicity of the 
control field is transferred to the logical Hamiltonian H{t), 
and an exact representation of the controlled evolution in 
terms of average Hamiltonian theory exists [6], [10], 



Ult) = 



7(k) 



iHt 



oo 

H^Y.-^'^' 



(6) 



each term H being computed from the Magnus series 
for H{t). As it turns out, the logical formulation is also 
particularly useful in situations where the control strategy 
directly incorporates symmetry criteria. 

For a realistic open quantum system, the influence of 
the surrounding environment may modify the dynamics 
in two important ways, (i) S may couple to a classical 
environment, effectively resulting into a (possibly random) 
time-dependent modification of the system parameters, in 
particular Hs ^ Hs{t). Deterministic time-dependent 
quantum control systems have been recently investigated 
in [15]. {ii) S may couple to a quantum environment E, that 
is a second quantum subsystem defined on a state space He 
of dimension >> ^5 and characterized by an internal 
Hamiltonian He- Let I5 ^ denote the identity operator on 
Hs,E, respectively. The drift Hamiltonian Ho{t) = HsE{t) 



of a general open quantum system may then be expressed 
as 

Ho{t) = Hs{t) (g,lE+Is(^HE+Y^ Ja{t) ® Ba , (7) 



where the Ba's are linearly independent environment oper- 
ators and, without loss of generality, we may assume the 
coupling operators (or error generators) to be traceless. 
In typical situations, both the exact time dependence of 
Hs{t) and Ja{t), as well as the exact form of He, Ba are 
unknown. If ps e (t) denotes the joint state of the composite 
S, E system, the evolution of S alone is now described by 
the reduced state obtained by a partial trace over E, 



Psit) = tTEipssit)) ■ 



(8) 



In general, the evolution of an initially pure state tts of S 
under the Hamiltonian Q, followed by (i) the ensemble 
average over the resulting time histories and/or (ii) the 
partial trace (|8}, results in a mixed state of S, tr(p|(t)) < 1. 
This implies genuinely non-unitary, irreversible dynamics 
for S, which physically accounts for quantum decoherence 
and dissipation effects [4]. 

For an open system, a control problem /orwa/fy similar 
to Q may still be formulated for the combined propagator 
U{t) of S plus E, provided that the action of the controller 
is explicitly restricted to the system variables only that is. 



H,{t) = H,{t)(E)lE, Ucit) = U,{t)®li 



(9) 



Two frame transformations may be relevant in the open 
system context. The transformation to a logical frame, 
which explicitly removes the applied control Hamiltonian, 
is effected as before, 

PSE{t) = Ul{t)psE{t)Uc{t) , (10) 

leading to a control problem formally similar to (0}, with 

HsE{t) = [U^,it)Hsit)U,{t)]®lE + Is®HE + 

+ [Ulit)Ja{t)Uc{t)] ^Ba . (11) 

a 

If a formulation which also removes the evolution due to 
He is needed, a simultaneous canonical transformation to 
a logical interaction frame is effected on the environment 
variables, 

p'sEit) ^ Ul{t)psEmEit) , UE{t)^e-^"-' . (12) 

The corresponding propagator U'{t) still satisfies an equa- 
tion similar to 0, where now 

H'sEit) - [ul{t)Hs{t)U,{t)]®lE+ (13) 

+ Y.Pl^^)'^-^'^)^M®[uUt)BaUE(t) . 

a 

The various propagators are related to each other as follows: 

U{t) = Uc{t)U{t) = Ucit)UEit)U'{t) . (14) 



B. Control tasks and performance indicators 

A dynamical control problem may be regarded as a steer- 
ing problem for the evolution operator of the target system 
in the appropriate frame. For an open system, a task of 
critical importance is decoherence control, which effectively 
requires the suppression of the error generators Ja{t). In 
particular, a decoupling problem consists in determining 
a control configuration {Hi,U(,{t)} such that for a given 
evolution time T > the joint propagator factorizes e.g.. 



U{T)^Xs{T)®Ue{T) 



(15) 



in the logical frame, Xs{T) being a unitary operator on 
S. Notice that Eq. M5\ implies decoupling in the physical 
frame as well. The simplest decoupling objective, on which 
I will focus henceforth, corresponds to identity design 
on S (the so-called no-op gate in quantum computation 
terminology [3], or complete decoupling or annihilation in 
decoupling terminology [6], [7]), whereby 



U{T) 



Ue{T) 



(16) 



If both Hs and the J a are constant in time, and Uc{t) is 
periodic, then the logical Hamiltonian (II It is also periodic 
and the above equation, once fulfilled at time T ~ Tc, 
remains valid for arbitrary times — NTc, N E "N. 
Under these conditions, the logical and physical frames 
overlap for every N, and the controlled evolution reads as 



Ps{Tn) = PsiTN) = Ps(0) ^TTs ^ Wii'l 



(17) 



Thus, arbitrary initial states of S are stroboscopically pre- 
served in both the logical and the physical frames. If either 
Hs or Ja are time-varying, and/or the control strategy is 
acyclic, it is still meaningful to require that 



PsiT) = ps{0) = TTs , T > 0, y-Ks 



(18) 



For stochastic control, the above objective is further relaxed 
to average state preservation in the logical frame that is, 

¥.{ps{T)} = ps{Q) = ns, T>0,y7Ts, (19) 

with E{ } denoting ensemble expectation. Clearly, con- 
trol schemes involving random operations are intrinsically 
acyclic, the control path practically never returning the 
system to the physical frame. If, however, the past control 
trajectory is recorded, this may be exploited to bring the 
state of S back to the physical frame at any time if desired. 

In order to quantify the accuracy of a given control 
procedure at achieving the intended objective, suitable per- 
formance indicators are needed. Let tt^ = I5 — 
denote the orthogonal complement of tts in Hs- Then the 
above task ( I19> is achieved if and only if, on average, the 
logical (reduced) state of the system has zero component 
along TTg (irrespective of the state of the environment). This 
naturally suggests to consider, for each pure initial state tts, 
the following a priori error probability. 



Note that eriirs) > for all tts follows from the fact that 
both TTg and ps{T) are Hermitian semi -positive definite 
operators. A worst-case pure state error probability may 
then be defined by maximizing over pure states that is. 



er = Max^sews {eTiirs)} 



(21) 



C. Control assumptions and group-theoretical design 

Control design is strongly influenced by the class of 
available controls. A particularly simple scenario is pro- 
vided by so-called quantum bang-bang controls [5], [6], 
whereby the control inputs Ui{t) are able to be turned 
on and off impulsively with unbounded strength, so as to 
implement sequences of effectively instantaneous control 
pulses. While such idealized assumptions must (and can [8]) 
be significantly weakened for realistic applications, the 
bang-bang setting provides the most convenient starting 
point for discussing stochastic schemes. 

Pictorially, it is helpful to visualize a control protocol 
in terms of the path that Uc{t) follows in the space of 
unitary transformations on S. For bang-bang controls, such 
a path is described as a piecewise constant time dependence, 
with jumps between consecutive values corresponding to the 
application of an instantaneous control kick. In particular, 
a large class of decoupling schemes may be obtained by 
constraining such values to belong to a discrete subgroup Q 
of unitary operators, the so-called decoupling group [6]. Let 
G = {gi}, where ge, £ ^ 0, . . . ,\g\ - 1, go = Is, denote 
group elements Cyclic decoupling according to Q over 
Tc is implemented by sequentially steering Uc{t) through 
each of the \Q\ group elements that is, 

Uc[{3-l)M + s\=g, , se[0,Ai), (22) 

, |^|. One can prove that. 



with M = Tj\g\ and j = 1 
in a fast control limit where 



Tc^Q, M ^00 , T = MTc > , 



(23) 



the leading contribution to the average Hamiltonian result- 
ing from HsE{t) in Eq. (II It is given by 

HsE = Hs (^Ie + Is He + ^Ja(g> Ba , 



X = ^ r dtUl{t)XUc{t) . 



(24) 



The advantage of group-based decoupling scheme is that 
the above time averages are directly mapped, via Eq. ( I22t . 
to averages over the control group Q, effectively implying 
a symmetrization of the controlled dynamics according to 
Q [6], [17], [16]. If, in particular, the action of Q is 
irreducible, then by Schur's lemma 



- 1 ^ tr(X) 



15 = 0, 



(25) 



eT(7rs) -E{trs 



' I am identifying an abstractly defined decoupling group with its image 
under a projective representation in Tis- Loosely speaking, 5 is a "group 
(20) up to phase factors", in general. This is irrelevant for the present discussion. 



immediately implying complete decoupling as in Eq. il6i . 

While cyclic schemes may be very powerful and concep- 
tually simple, they are only applicable (at least in the simple 
formulation presented here) to time-independent control 
systems. Also, because averaging requires traversing all of 
Q, they tend to become very inefficient as the size of Q 
grows. The basic idea that underlies random decoupling 
according to Q is to replace sequential cycling with random 
sampling over Q. In the simplest kind of protocols, the value 
of the propagator Uc{t) is determined by a group element 
which is picked uniformly at random in Q that is. 



Prob {gi 



\Q\ 



(26) 



Thus, both the past control operations and the times at 
which they are effected are known, but the future control 
path is random. Under these conditions, no average Hamil- 
tonian formulation is viable, and averaging effects emerge 
through ensemble rather than time averages. 



{{X{t)))^^{ul{t)X{t)U,{t)] 



(27) 



Under the uniformity assumption, such expectation values 
again reduce to averages over Q, leading to the possibility 
of stochastic averaging. 



(28) 



The two key questions to address for random decoupling 
are to understand whether stochastic protocols are indeed 
capable of achieving decoupling and, if so, how they per- 
form compared to deterministic counterparts. We focus here 
on the first question, by presenting an explicit derivation of 
an error bound for randomized control directly within the 
open-system context ^. 

III. Random decoupling 

A. General error bounds 

We begin by recalling a few preliminary facts. 
Remark 3.1: Let \\A\\2 = Max |eig (a/aTA)] denote the 
operator 2-norm of A. Then (see e.g. [18]) 

ii) ||A||2 =Max|eig(A)| , VA = At; 

\\AB\\2<\\AM\B\\,, yA,B- 
(m) If U is unitary, \\WAU\\2 = \ \A\\2 , VA 
Lemma 3.2: Let A be any rank-1 operator on Tis- Then 



MA)\<\\A\\2 

Proof: A may be represented as A 
dimensional complex vector \v) = [vi,. 
lbll=Efekf-Then 



|ti)(w|, for a d- 
, Vd] with norm 



|tr(A)| = \v,\ < ||«|| =Max|eig(/^)| = \\A\\2 . 
Q.E.D. 

~In [13], a detailed proof was obtained for the closed-system setting, 
and used to sketch the main steps leading to the open-system result. 



Theorem 3.3: Let S be an open quantum system de- 
scribed by a Hamiltonian of the form 0. Suppose that the 
control protocol satisfies the following assumptions: 
(j) (Irreducibility) Q acts irreducibly on Tis- 
(a) (Uniformity) Uc{t) is uniformly random for each t. 
(Hi) (Independence) For any t,s > 0, Uc{t) and Uc{t + s) 
are independent for s > At. 

If, in addition, the total interaction Hamiltonian is uni- 
formly bounded in time. 



Hs{t) 



Ja{t) (g) Bait) 



<k, yt, (29) 



then 



€T = (TAt fc2) for TAt < 1 



(30) 



Proof: Let tts be an arbitrary pure state of S. The first 
step is to cast the pure-state error probability (I20> in a more 
convenient form to bound. By purifying the initial state of 
E if necessary, we may assume that pse{0) — tts ® t^e, 
both TTs,E being one-dimensional projectors. By using the 
definition of partial trace and the cyclicity property of the 
full trace, we have 

eriiTs) = E{trs(4ps(T))} 

= E {tisE (tts 

= E jtrsB (tt 

= E jtrsB (tt 

where the relation (I14> has been used, and UE{t) drops. 
Let H'g^{t) denote the interaction Hamiltonian of Eq. M9\ . 
Then the task is to bound the error in implementing identity 
design on the logical interaction propagator at time T, 



4' 



(31) 

Ie~Pse{T))] 

1eU{T)tts (Xi tte 
lEU'{T)ns'E>7TEU'{T)^)} , 



U'{T) = Texp <^ / duH'sEiu) 



(32) 



with H'sEit) ^ Ul{t)H'sE{t)Uc{t) given in Eq. O- 
The above propagator may be expressed as follows: 



(7'(T) = ^/„(T), 



(33) 



n=0 



/„(T) = (-Z)" / dMH'sE{Un)...H'sE{ui). 
Jo<ui...<u„<T 

(34) 

and similarly for U{T)\ with du = dui . . . dun. Thus, we 
need to calculate 

(tts) = 

= E|tr5£( ^7rs®^i5/™(T)t4®lB/„(r)jl . 

I \n,m=0 /J 

The contributions with n = or m = vanish because of 
TTg and TTs cancel each other upon exploiting the cyclicity 
of the trace. Because triT^s) > 0, 

\^t{t^s)\ < 

< J2 |E{trsB(^5 ® 7rB/„(r)t4 ® Ii;/„(T)) | 



Under the assumption of sufficiently smooth behavior, the 
expectation may be moved under the integral. Fix a pair of 
integers n,m > 1, then the relevant contribution is 



dudt E (tTs ® TTEH'sEitl) ■ ■ ■ H'sE{tm)Trs ® 



H'sE{un)...H'sE{ui)] ,(35) 



where the integration region VK^"^™) — {(u, t)|0 < ui 
. . . < u„ < T; < . . . < i„ < T}. Let w[''''^\l^t) C 
I^C".™) denote the subset of points satisfying that ui, tg 
are each time-ordered and no ui or tg is further away than 
At from the rest, and let W^'''"'\M) C denote 
the remaining region. Because, within Wj"'™'' (Ai), at least 
one of the integrating variables is more than At away from 
all the other variables, the independence assumption {in) 
allows the expectation relative to such a variable to be taken 
separately. By the uniformity assumption (m) on Uc{t) for 
all t, and by the tracelessness assumption on HgEit) for all 
t, such an expectation vanishes. Therefore, W'f'''"'' (At) is 
the only subset of points contributing to the expectation in 
Eq. (I35> . Let dw^"'™^ denote the corresponding integration 
measure. Then 



E 



n,m>l 1 

{trs£;( 



'se{ui))] 



< 



E 



W- 



E. 



{|trsB( 



n.'m> 1 ' 



IE ■ ■ 



.h: 



seM)\} 



where in the second step Jensen's inequality has been used. 
By noticing that the argument of the trace is a rank-1 
operator. Lemma I3T2I mav be used to simplify 



E 



eriTTs) < E / 

{ 



n,m> 1 ' 



TTS (g) TTEH'sEih) . ■ - TTg <^Ie ■ ■ ■ -^SE: ("l ) | | ^ } 

< E Vol(VFi("''"^)fc"+" , 



where the inequality [ii) in the Remark 13.11 and the 
uniform bound k for H'g^{t) in M9\ have been used, 
and Vo^M^i^"'™^) is the volume of w[''''^\ Note that the 
dependence upon Trg has disappeared at this point. 

The above volume may be estimated through a combi- 
natorial argument. First, notice that given the two ordered 

lists Q < Ui < . . . < Un < T, Q < ti < . . . < Un < T, 

there are ("^™) different merged orderings. Fix a particular 
one. Then each element needs to be either within At of the 
next one or of the previous one. Make a choice for the 
odd-numbered elements, the first element being labeled 1. 
There are at most 2r("+™)/2T such choices. For each of 
them the contribution to the volume may be bounded by 



ordering the even-numbered elements, then by inserting the 
odd ones, ignoring the ordering constraint now. Finally, 

^1 ' - \ m ) ([(n + m)/2j)! 

2\(n+m)/2\j,l(n+rn)/2\ r("+™)/2l 

= Km , (36) 

where the inequalities ("+") < 2"+"-i (for n + m > 2), 
and [{n + m)/2j! > 2L("+'")/2J-i have been exploited. 
The last step is to sum over n, m: 



Max,,{eT(^5)} = eT < E 



(37) 



This may be done by considering separately the four partial 
sums where both n and m have the same (even or odd) 
parity, or they have opposite (even-odd or odd-even) parity, 
respectively, and by evaluating the [J, [] in Eq. i36i 
accordingly. Lengthy but straightforward calculations yield 

..<(4TAtfc^^l + «^^^ + 4^^^^'-— '^^ 



0{TAtk^) , 

(38) 



(1 - ATAtk-^Y 



for values of TAtk^ ^ 1, as quoted in Theorem 13. 31 
Q.E.D. 

Remark 3.4: By setting all the coupling operators Ja = 
0, the error bound for random decoupling of a closed or 
classically time-dependent control system is obtained. 

According to the above Theorem, the performance of 
stochastic control can be made arbitrarily high by appro- 
priate design, in particular by choosing a sufficiently small 
At in the present setting. Remarkably, this implies the 
possibility to arbitrarily suppress on average decoherence 
in the logical frame. Note that, unlike deterministic decou- 
pling, stochastic schemes place no restriction on the time 
dependence of Ho{t), only on the maximum eigenvalue 
of the interaction part, Hg^{t). The latter, however, may 
diverge in physical situations involving infinite-dimensional 
environments. Thus, appropriate care is needed to properly 
define the relevant strength k in such situations [19], [20]. 
Physically, the parameter k~^ is of the order of the shortest 
correlation time present in the interaction to be removed. 
While this provides the relevant time scale to the purposes 
of obtaining an upper error bound, lower or typical error 
bounds may be better in specific situations, depending on 
the details of both the system and the environment. 

B. Example: Control of a single noisy qubit 

A simple illustrative example is provided by a single two- 
state system (a qubit) dissipatively coupled to a quantum 
reservoir. In this case 71$ — and a basis for the traceless 
operators on S is given by the Pauli operators, a^, a = 
x,y,z. Consider for simplicity a time-independent open- 
system dynamics. Eq. takes then the form 



Ho = LjQ(7^ ®Ie + ls ® He + y^^(Jo 



(39) 



where cr^ represents the energy eigenbasis of the isolated 
qubit, and oja, Ba are appropriate real parameters and 
Hermitian environment operators, respectively. Complete 
decoupling may be achieved in the deterministic setting 
by cycling the control propagator through a (projectively 
represented) error group for the qubit^* that is, Qp ~ 
{Is, <Jx, (Jy, cFz}- Thus, Tc ~ 4At, and Eq. (I22> yields 

f Is ieAti, 

[ t <E Ati . 

In practice, this corresponds to a series of four equally 
spaced bang-bang so-called tt- (or 180°-) pulses, alternating 
between the x and z axes. In terms of the control inputs 
introduced in Q, a vr-pulse along the a axis may be 
performed by applying a linearly polarized oscillating field 

HaUa{t) = CTq y{t) COs[ljj(t — tp)] , 

v{t) ^v[e{t~tp) -e{t~tp -t)] , v>o, 

where oj — ujq on resonance, tp, t aie the time at which the 
pulse is applied and its duration, respectively, and 2Vt — n 
with T — > 0, y ^ oo to satisfy the bang-bang requirement. 

For random decoupling over the Pauli group Qp, the 
control prescription i26\ corresponds to applying a sequence 
of TT-pulses with are randomly drawn from Q p that is, each 
of the Pauli operators is applied with probability 0.25 at 
times tj = jAt, j S N. Physically, the relevant strength 
parameter k may be associated to the high-frequency cut- 
off uJc that is contained in the reservoir power spectrum 
and determines its frequency response. In general, however, 
additional time scales related to both uiq and the temperature 
affect the overall control performance. Thus, according to 
the worst-case bound of Eq. ( I30t . decoherence suppression 
at time T is achieved provided At is made sufficiently small 
with respect to uj^^. Remarkably, an exact solution for the 
stochastically controlled dynamics may be obtained in the 
special case where = By = 0, corresponding to pure 
decoherence. A detailed analysis of this limiting situation 
is reported in [21]. 

IV. Conclusion 

I have discussed a control-theoretic formulation which 
explicitly invokes random control design, and which is 
applicable to arbitrary finite-dimensional, time-dependent 
open quantum control systems. I focused on random de- 
coupling design for decoherence suppression as a relevant 
case study, and showed how arbitrarily low error rates may 
be achieved in principle. Further study is needed to both 
explore concrete applications of randomized schemes and 
assess their full potential, as well as to integrate random 
design within existing control settings. Beside pointing to a 
still largely unexplored territory in the theory and practice 
of quantum control, the ideas presented here might allow to 

^^The abstract decoupling group is Z2 X Z2 in this case. 



take advantage of novel perspectives, as offered for instance 
by noisy quantum games [22] or randomized algorithms 
for classical uncertain systems [23]. It is my hope that 
the results presented here will prompt the control theory 
community to further investigate the interplay between 
randomness and coherence in quantum dynamical systems. 
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